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Parity violating electron-deuteron scattering can potentially provide a clean access to electroweak 
couplings that are sensitive to physics beyond the Standard Model. However hadronic effects can 
contaminate their extraction from high-precision measurements. Power-suppressed contributions 
are one of the main sources of uncertainties along with charge-symmetry violating effects in leading- 
twist parton densities. In this work we calculate the twist-four correlation functions contributing to 
the left-right polarization asymmetry making use of nucleon multiparton light-cone wave functions. 
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I. INTRODUCTION 

Even after decades of experimental studies, deep in- 
elastic scattering (DIS) remains one of the most power- 
ful tools for unraveling the partonic structure of nucleons 
and nuclei. DIS also allows for systematic searches for 
physics beyond the standard model. Parity violation in 
DIS (PV-DIS) at medium energies is particularly sensi- 
tive to effects of New Physics. Historically, this process 
played an important role in verifying the Standard Model 
d Q. Today the search for New Physics motivates a 
number of ongoing and planned experiments [3- 9]. The 
physics reason for this great interest is that within the 
standard model the Weinberg angle Ow should show a 
highly non-trivial characteristic scale dependence, which 
can be mapped out by combining experiments at differ- 
ent momentum scales. The SoLID experiment at JLab 
[Tol fllj (see also [IH, [l]|) will be especially sensitive to 
the poorly measured weak neutral coupling constants Ci q 
in the low-energy electroweak Lagrangian 
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To analyze the theoretical situation, effects of New 
Physics are parameterized by SCi a [46| according to 
Ci a = 2g A g v +SCi a and C 2ct = 2gyg%+SC2 a , where the 
standard model coupling constants are g v A = Q^^Q^j 
in terms of the left and right (a = L, R) weak charges 



QZj=T 3 (f a )-Q(f) sin 2 9 
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The (SCi a ) are inaccessible in other measurements, which 
gives PV-DIS its unique quality. 
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The projected sensitivity of the SoLID experiment for 
an asymmetry discussed below is SA/A — ±0.005(stat.) 
at an average Q 2 of 3.3 GeV 2 and an average x of 
(a;) = 0.34, which sets the scale for the size of accept- 
able theoretical uncertainties. At this level of precision 
several sources of systematic uncertainties can hamper a 
precise determination of the Ci a , as discussed recently in 
Refs. ITU . Some of the most relevant are uncertainties 
in leading-twist parton distributions functions, in par- 
ticular charge-symmetry violation (CSV), contributions 
from higher- twist correlation functions, and kinemati- 
cal target-mass corrections. Far from being a nuisance 
higher-twist correlations encode very interesting and yet 
little known information on hadron structure. Therefore, 
all cases in which leading-twist contributions are absent 
or reduced, such that one has a good chance to determine 
higher- twist ones are of great interest. If the relevant 
higher-twist contributions are measurable with a given 
experimental sensitivity, as we will claim they are not in 
this case, one is in a win-win situation: PV-DIS can be 
regarded either as a tool to find New Physics, in case the 
effects of the latter are prominent, or it can be seen as 
a venue to access unknown aspects of strong interaction 
physics. 

Parity violating weak interactions give rise to an asym- 
metry in the inclusive cross sections for scattering of left- 
and right-handed electrons off a deuteron 
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This is the main medium-energy observable in PV- 
DIS which will be scrutinized at Jefferson Lab [Tol - FT^ j . 
Among all uncertainties of the theoretical prediction of 
this asymmetry we will focus on the power suppressed 
contributions. Two recent studies of it reached some- 
what different conclusions (T3 . Il5|. Our results turn out 
to be very similar to those from [14|. 

It was demonstrated by Bjorken and Wolfcnstein 
[l6L that twist-four corrections to the asymmetry 
are due to a single (nonlocal) four-quark operator. The 
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FIG. 1: Kinematics in deep inelastic deuteron electron scat- 
tering. 



first estimates of the matrix element of the spin-two part 
of this operator were obtained in the framework of the 
MIT bag model flU [Hj]. This technique was extended 
in Ref. [14| to include the effects of higher spin opera- 
tors. It was found that their effect is negligible within 
the model used. Renormalon analysis offers yet another 
technique to model the momentum fraction dependence 
of certain higher twist matrix elements |20L |21| . These 
renormalon-based studies demonstrate [2(| that higher- 
twist correlation functions (involving two quarks and a 
gluon) tend to grow at large x, i.e., like (1 — x)~ , in 
qualitative agreement with experimental measurements 
of electroweak structure functions [22|. However, the 
four-quark operators we consider are free from ultravi- 
olet renormalons [23| and thus this approach is not ap- 
plicable. The absence of renormalon contributions might 
explain the qualitatively different behavior of such corre- 
lators and gluonic ones. In this work we calculate twist- 
four corrections employing a model for the nucleon wave 
functions in the light-cone formalism which was proposed 
by [Hl|. 

The paper is organized as follows: Sect. HI] contains 
basic definitions and notations. In Sect. IPlTl we give a de- 
tailed discussion of power corrections to the asymmetry 
(|3]). In Sect.|TV]the necessary ingredients of the light-cone 
formalism are given. Results of our calculation and our 
prediction for the twist-four corrections to the asymme- 
try are collected in Sect. [V] Finally we give our conclu- 
sions. Several Appendices contain technical details and 
formulae left out in the body of the paper. 



II. PRELIMINARIES 

Let us briefly discuss the physical observables we will 
be analyzing below. The cross section for polarized elec- 



tron scattering off an unpolarized deuteron target, with 
kinematics shown in Fig. [TJ is given by the sum of three 
terms 



da L / R = do L J R + dotl* + 2do h J* 



(4) 



which describe the contributions due to the electromag- 
netic and weak interactions and their interference. Each 
term is a function of the standard kinematical variables 



-q 2 , v = (p-q), x 



2v 



2 »=M, (5) 



Each term in Eq. Q is given by the convolution of a lep- 
tonic and hadronic tensor. This reads in the laboratory 
frame 



da 



L/R 



where the repeated Latin indices imply summation over 
electromagnetic and weak exchanges a,b = (e,w). The 
coefficients 
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encode the products of gauge boson propagators and in- 
teraction strengths. The leptonic tensor admits the con- 
ventional form 



(T L/R\ _ n L/R n L/R p L/R 
\^ab IV-v — Vo "-av l 



(7) 



which is a product of the electromagnetic (weak) charge 
S3 Qe(w) f° r tne ^ (right) handed electron 



Q 



and 



L/R 



Qw = -o +sm 2 e w , Q* = sm 2 6 w , 



t v — kpk v + k^k v g^ v (k ■ k ) i is ^vpak^k . (8) 

The hadronic tensor is the deuteron matrix clement 
of the product of currents 

1 



d 4 ze iq - z 



8ttM d 

x (D(p)\{j£(z)%(0) +%(z)%{0)}\D(p)) , (9) 

where Mu is the deuteron mass and averaging over 
deuteron polarizations is implied. The electromagnetic 
and neutral quark current are defined as (cf. Eq. (|2"]l) 

f e = qQfq , & = q L r 3 rqL - sin 2 6 w f e , (10) 

where q — (u,d). It was demonstrated by Bjorken (l6j 
that if one assumes valence quark dominance in the re- 
gion x > 0.4 and neglects all sea quark and isospin break- 
ing effects (which should be justified for large virtual 
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mass Q 2 ), the asymmetry ([3]) becomes free of hadronic 
physics contaminations and is given by the Cahn-Gilman 
formula 



the form 
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New Physics is best parameterized by allowing for non- 
standard values for the coefficients Ci a , which are rein- 
troduced in Eq.ljllj) by replacing 
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20 _• 2 n 
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1-4 sin 2 9 W -> 



(2Ci u — C\d) , 

(2C2« — C2d) • 



However, the assumptions leading to vanishing 
hadronic effects are only valid approximately and have 
to be abandoned in the analysis of high precision experi- 
ments. The main hadronic effects are caused by CSV and 
power suppressed correlators. The central point behind 
our work, and that of others, is that these effects have a 
strong x dependence which allows, if precisely known, to 
isolate and subtract them and thus to increase the sen- 
sitivity of experiments like SoLID to New Physics. Thus 
one has to go beyond leading approximations and has to 
take into account higher-order electro-weak effects, sea 
quark effects, target mass and higher-twist corrections 
at least at a level matching the accuracy of experimental 
measurements. 

CSV arises from isospin violation of u and d quark 
distributions in the proton and neutron, i.e., by Su = 
Up u n ^ and 5d = d p — u n ^ 0. Modern global anal- 
ysis of parton distribution functions incorporate CSV ef- 
fects, which are found to become more significant as x 
decreases [13], R C sv ~ (Su — Sd)/(u + d) ~ (1 — x) 4 ^/x. 
CSV effects might explain a significant fraction of the dis- 
crepancy between the NuTeV results [28| and predictions 
based on the standard model and isospin symmetry. 

The other source of corrections are power suppressed 
contributions from multi- particle correlation functions. 
Obviously the nucleon wave function is a complex state 
containing many highly entangled Fock states, only par- 
tially characterized by parton distribution functions. The 
isolation and determination of specific multiple-field cor- 
relators is the logical next step to explore hadrons and 
is therefore of great interest in its own right. In con- 
trast to mere one particle probability distributions, they 
contain information on relative phases. As they are typ- 
ically power suppressed, high luminosity experiments at 
medium large Q 2 are needed to extract them. These are 
requirements which are perfectly fit by Jefferson Lab, es- 
pecially after the energy upgrade. 



III. TWIST FOUR CORRECTIONS 



In the region of low transferred momentum Q 2 <C M 2 ^ 
one has da ww <C da ew <C da ee and the asymmetry takes 
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where we took into account that da^ e = da^ e = da ee . 

Introducing the scalar, isovector and axial isovector 
currents 



q^T 3 q , 



A" 



qri5T 3 q (13) 



one can represent the electromagnetic (weak) hadronic 
tensors as follows 

W£ (p, q) =WT (P, q) + ~W% v [p, q) , 



- sin 2 e)WP(p,q) 



Ism 2 eWg v (p,q)-lw^(p,q), (14) 



where 



W^{p,q) 



W^(p,q) 
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d 4 ze lq - z (D(p)\V f "(z)V u (0)\D{p)} , 

d 4 ze^ z (D(p)\S»(z)S»(0)\D(p)), 

d 4 ze iq - z (D(p)\A"(z)V v {0) 

+ V»(z)A»(p)\D(p)}. (15) 



Here we took into account that the deuteron matrix ele- 
ments of nonsinglet terms, i.e., involving the product of 
isovector and isosinglet currents VS and AS, vanish by 
isospin symmetry since the deuteron is an isoscalar state. 
Keeping only twist-two terms in the OPE expansion of 
the hadronic tensors (|T5|) one arrives at the Cahn-Gilman 
formula (llip . the first and the second term in the square 
brackets in (fTTj) arise from vector- vector (W^ w ) and axial- 
vector (Wew) correlators, respectively. The corrections 
to the Cahn-Gilman formula can be parameterized as 
follows 



A RL = - 



G F Q 2 3 



ai + a 2 



i-(i-y) 2 

1 + (1 - y) 2 



(16) 



where (i = 1, 2) 



= -(2C lu -C ld )[l + R l 



(17) 

Here the functions Ri (i = 1, 2) alluded to before receive 
contributions from several sources of hadronic effects. 
The precision measurement of the mixing angle at low Q 2 
gives sin 2 9w — 0.2397 [5|. Thus the axial current con- 
tribution (cb) to the asymmetry is relatively small and 
we will focus on the calculation of twist-four corrections 
to 5i. They can be easily identified. Indeed, neglecting 
effects of isospin breaking one gets (see Ref. |la |) 

(D\S M (z)S v (fl) ~ V»(z)V v (0)\D) = 

= ^(D\u{z)^u{z) d(0)Yd(0) + (u^ d)\D) . (18) 
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The expansion of the operator at the right-hand side of 
this equation starts from twist-four. In terms of 
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d 4 ze iq - z 



(19) 
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x (D(p)\u{z)^u{z) d(0)j v d(0) + (u d)\D(p)) 

we define the structure functions F^ =1 2 as coefficients in 
the tensor decomposition 

M D wr(p,q) = (-g^- 



q^q v 



pa 



(20) 

where the index runs over a = V, S, ud. The twist-four 
contribution to R\ takes the form 

pud 



1 



10(1 
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where 



T a = xyF? - 



1 - 



, sin 2 (V) T s 
1 xM D 



u 



2E 



F% 



(21) 



(22) 



Keeping in J- s and J- ud the dominant contributions only, 
i.e, twist-two and twist-four, respectively, and taking into 
account that they both satisfy the Callan-Gross relation 
F 2 = 2xFi , one finds 

pud pud 

WW' < 23 » 

The expression for Ff at lowest order of perturbation 
theory is given by the sum of parton densities in the 
deuteron 

F i ( x ) = I \ u d ( x ) + d D (x) + u D (x) + d D (x)] , (24) 

where as usual q~D(x) = — Qd(~ The quark distri- 
bution functions are defined by the matrix elements of 
nonlocal light-cone operators, 

(D\q(z)itq(-z)\D) = 2(p-z) [ dxe 2l ^ x q D {x) . (25) 



To evaluate F^ d we represent the hadronic tensor W£j 
via the dispersion relation as a time-ordered product of 
electrowcak currents 



d 4 ze iq - z 



(26) 
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x (D(p)\T{u(z)^u(z) d(0)Yd(0) + (u^ d)}\D(p)) 
and make use of the operator product expansion [29| 

r _ >, tw — 4 

T< u(z)jfj,u(z) d(— z)j u d(— z) + (u f-> d) > 



, logz d^dv j du—Q(uz) 

+ \s^ z a d? f —Q(uz)\, (27) 
z Jo u ) 



where S uav p = gua9vfi + g va g^i3 - 9u V g a p- 

The operator Q (Q 2 m the notations of Ref. [29|]) is 
given by the following expression 



Q(z) =i J dv J dt 



nr 2 n 3 - 4 QWM,i,-i) 



n+n+ Q A (i, v ,t,-i) 



+ (z <-> -z) . (28) 



Here 



Qa{o) =(u(aiz)t a ^ 5 u(a 2 z)j (d(a 3 z)t a ^d(a±z) 
Qv(o) ^(u(a 1 z)t a ^u(a 2 z)^ (d(a 3 z)t a ^(a^zj) , (29) 

and Tlf k = (1 ± Pik), where Pik is the permutation oper- 
ator, e.g., Pi 2 Qy(ai,a 2 ,a 3 ,a 4 ) = Qv (a 2 , ai, a 3 , a 4 ). For 
later convenience we rewrite (1251 as follows 



Q(z) =i j dv J dt 



e+(i,w,t,-i)-c-(i. «.*.-!)]. 

(30) 



where 



Q+{a) =(1 + Pi 2 P 34 )(l + Pi 4 P 23 )Q+(a) , 
C_(o) = (Pi 2 +P 34 )(l + Pi 4 P 23 )Q-(a) (31) 



and 



Q±(a) = Q v (a)±Q A (a). 



(32) 



Let us define the twist-four distribution Qd{x) as a 
deuteron matrix element of the operator Q 

(D\Q{z)\D) = i [ dxe 2i ^ x Q D (x). (33) 



It follows from ((28)) and (|32|) that Qi)(a;) is an even func- 
tion of x with vanishing first moment, 



dx Q D (x) = 0. 



Inserting (|27|) and ([28|) into (|26l) one finds after some al- 
gebra 



Ff{x) = -^xQ D {x). 



(34) 



Then, keeping in Ff the valence quark contribution only 
we obtain the following expression for the twist-four cor- 
rection to the asymmetry 



1? 



x Qd(x) 



Q 2 5(1 - f sin 2 Qw) u D (x) + d D {x) 



(35) 



The deuteron is a weakly coupled state of the proton 
and neutron with the binding energy Eb — 2.2 MeV. In 
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the incoherent impulse approximation its hadronic tensor 
in the deuteron's rest frame can be represented as [30j 



W D (p a) ~ / ^ 



\f(Ps 



(36) 



where p s = (E Pa ,p s ) and the integration is performed 
over the spectator three- momentum p s , see Fig. [TJ Here 
f(jp s ) is the deuteron wave function in its rest frame, 



normalized as 



1 



and Whu 



(p(«)) 



(27r)-3fd 3 p s \f(p s )\ 

are the proton (neutron) hadronic tensors. The function 
f(p s ) is strongly peaked at p s = [30]. Thus one can 
simplify the above expression by neglecting terms of order 
~ \p s \/Mjf and higher under the integral. Then one finds 

W£ (p, q) c W$ (p/2, q) + W$> (p/2, q) , (37) 



and as a consequence G?Od ~ do v + c?cr„. Then Eq. (|37|) 
yields the following relation between the structure func- 
tions of deuteron and nucleons, 

Fi(x/2)^F%(x)+F?(x). 

It turns out that this approximation overestimates the 
deuteron structure function by 5 4- 10% [13, [HJ. This is 
acceptable for our purposes. For the parton densities the 
corresponding relation reads (cf. Eqs. (47) and (48) in 
Ref. O) 



-q D (x/2) ~q p {x) + q n (x) . 



(38) 



Similarly, defining the proton (neutron) twist-four distri- 
butions by 



(N\Q(z)\N) = i / dxe 2i ^ x Q p {x) 



(39) 



one gets for the deuteron twist- four function Qd{x) 

\Qd(x/2) =Q p (x) + Q n (x) = 2Q p (x) . (40) 



IV. NUCLEON LIGHT-CONE WAVE 
FUNCTIONS 



Our lack of information on the magnitude of higher- 
twist matrix elements is the main obstacle for a quan- 
titative analysis of power-suppressed contributions to 
hadronic cross sections. Hadron structure models provide 
estimates for the size of nonperturbative matrix elements, 
but their predictions vary strongly. This is understand- 
able in view of the fact that confinement is incorporated 
rather differently The first estimates of twist-four cor- 
rections to the asymmetry ([3]) were obtained within the 
MIT bag model [lg, [l9( which incorporates confinement 
quite ad hoc, (see also Refs. [H], El] for recent devel- 
opments). In this work we use another approach, the 
light-cone formalism [Hj] , for the evaluation of twist-four 
corrections. 

In the light-cone formalism the nucleon is represented 
by a superposition of multi-parton Fock state wave func- 
tions. The latter are functions of the parton longitu- 
dinal momentum fractions Xi, transverse momenta fe_u, 
and parton helicities. The light-cone wave functions 
(LCWFs) are eigenfunctions of the QCD Hamiltonian 
quantized in the light-cone gauge [13, IHj]. Models for 
LCWFs of various degree of sophistication have been con- 
sidered in different context in the vast literature on the 
subject, see, e.g., Refs. [^.[2E.[3l[36r|38t. In this work we 
will follow the formalism developed in Refs. [24Tj26l l36j 
and will take into account only the lowest components of 
the nucleon LCWFs: the three quark and three-quark- 
gluon component. The details of the light-cone formal- 
ism relevant for our further discussion are collected in 
Appendix [X] 

The three quark component of the nucleon state is pa- 
rameterized in terms of corresponding LCWF as 
follows 



Here we took into account that Q P {x) — Q n {x) due to 
isospin symmetry. 

We also define the nucleon twist-four distribution 
Q±(£) (and similarly Q±(£)) by 

(N\Q±(a)\N) =(p-zf J 2?£e~' i(p - 2)Efcafc4fc Q±(0, 

(41) 

where £ cumulatively denotes the array of four variables 
£ = (£1, £2, £3, £4) and the integration measure stands for 

v £ = nLi d & S(Ei &)• Then it; follows from Eq. (|30]l 
that 



Q P (x) 



P£ 



-{(6 + + ft + 6) 



£ 2 £3(£2+£ 3 ) 

£ 3 <5(z + £1) - £ 2 <5(£ 4 - a;)}[Q + (0 - Q-(£)] • (42) 
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[DX] 3 ^ 3 (X) 



(4 T (l)ut 4 (2)4 f (3) - uf t (l)(4(2)<4 t (3)) |0) . (43) 



Here and below for notational simplicity arguments like 
I in u\Al), stand for the collection of all relevant argu- 
ments, i.e., itj,j(<?) = u|j(x£, k±e). The creation (annihila- 
tion) operators of a quark with helicity A and momentum 
p satisfy the commutation relation (IA.10[) . As usual, the 
momentum fraction Xi is defined as ratio of the longi- 
tudinal (i.e., "+") momentum of the i— th parton and 
the one of the nucleon. The integration measure has the 



6 



following form 



namely, 



[VX] N =- 



1 



\J%\ ■■■XN 
N 

[dx] N = dxi 6(1 - y^Xj 



dx]N[d 2 kA_] N , 



[d 2 k ± ] N 



(16tt 3 ) 



i N 

lyr^U/k^ ($>„ 



(44) 



Here we accept the Bolz-Kroll ansatz [24j for the func- 



tion * 



(0) 
123 



(45) 



The transverse momentum dependence is encoded in the 
function fijv 



Q N (ciN,Xi, fcj 



2^2 \JV-1 



(16^) 

X\X 2 ■ ■ ■ Xn 



■ exp 



(46) 



which is normalized such that 

[d 2 k±\ N £l N (a N ,Xi,k±i) = 1, 

[d 2 k ± ] N fl 2 N (a N ,Xi,k ±i ) = — 
J xi 



Pn 



..x N 



(47) 



where p^ = (87r 2 a^,) JV_1 . The function (j>(xi), entering 
(|45[) . depends only on the longitudinal momentum frac- 
tions of constituent partons and is related to the leading- 
twist, i.e., twist-three, nucleon distribution amplitude, 



J 



(xi,x 2 ,x 3 ) = $ 3 (xi,X2,x 3 ;fj,a) . 



(48) 



Here $3(2;) is the twist-three nucleon distribution ampli- 
tude defined at the low-energy scale po = lGeV. We use 
the following ansatz for $3(0;) [24[ 



$ 3 (x) = 60xix 2 x 3 (1 + 3xi) 



(49) 



which emerges from the truncation of the conformal par- 
tial wave expansion after the lowest few terms. The nor- 
malization constant /jy in Eq. (|45|) is determined by the 
matrix element of the corresponding local three-quark 
operator. The analysis within the framework of QCD 
sum rules [3{| yields in the following estimate for fjy 
[UEI at the scale p = 1 GeV 



f N = (5.0 ±0.5) x 10" 3 GeV 2 



(50) 



On the other hand, the parameter a 3 determines the 
smearing of the wave function in the transverse plane 
and, e.g., the average quark transverse momentum. Fol- 
lowing Ref. we take 03 = 0.73 GeV" 1 in our esti- 
mates. With this set of parameters, the contribution of 
the three-quark Fock state to the norm of the nucleon 
state is about 17%, 



P. 



3q 



435 
112 



ftps" 0.17. 



(51) 



The four-parton quark-gluon contributions with zero an- 
gular momentum to the nucleon states have the following 
form [H 



|.Pj ~\~)uudg± 



\P,+) 



uudg^ 



Jjk 



-VX] 4 *[ 234 (X)af(4) [t°n t (l)]tn J t t (2))4 t (3)|0) , 
-PXU{*IW{X) [*X(l)]f (4.(2)4^(3) - 4 t (2)4 4 (3))af(4) 

+ ^(X)ul(l){[t a H (2)}] 4 f (3) - [i%(2)]t 4 t (3))«? !t (4)||0) , 

I 



(52) 



where the four-parton LCWFs are again taken in the 



Bolz-Kroll form 
1 



^1234 



/2x4 



4> g (xi,X2,x 3 ^ xa) fli(a g ,Xi,k±i) , 



*1234 = m — ^ ^2, x 3 , Xi) Q 4 (%, Xi, kj_i) , 

\J iXi 

V&12I4 = -7= — ^1 2) (xi,X2, x 3 , xa) rLt(a g , Xi, k ±i ) . (53) 

J IXa 
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The functions ipg^ipg 1 ^ which depend on the light-cone 
momentum fractions of the partons can be expressed in 
terms of the twist-four quark-gluon nucleon distribution 



amplitudes introduced in Ref. 45| . 

g<pg(xi,x 3 ,x 2 ,X4) = 

2Hf (xi, X 2 , X 3 , Xi) + EI(X2,X1,X 3 ,X4,) 



- mN 
96 

gtfj g 1) (x 1 ,X2,X3,X 4 ) = 



m N 
' 48 



^\{xi,X\,Xl,,X±) + -<!>l(xi,X2,X 3 ,X4,) 



gip q 2) (xi,x 3 ,x 2 ,X4) 



m N 
48 



$l(xi,X 2 ,X 3 ,X4) + -^l{x 2 ,Xi,X 3 ,X4) .(54) 



Keeping only the lowest terms in the conformal expansion 
of the corresponding distribution amplitudes one arrives 
at the following expressions [26J 

g(/> g (x 1 ,x 2 ,x 3 ,x 4 ) = - 210m N X 9 1 x 1 x 2 x 3 xl , 

g% 1 '(x 1 ,X2,x 3 ,x i ) = - 105mjv(Af + A§) xix 2 x 3 x\ , 



g^f ) {xi,X2,x 3l x i ) = - 105TOAr(Af - Xl)xix 2 x 3 3 



(55) 



The sum rule technique was found to give the following 
estimates for the coupling constants X 9 k at low energy 
scale 1 GeV [H 

Af =(2.6 ±1.2) ■ l(T 3 GeV 2 , 
Af =(2.3 ±0.7) • !CT 3 GeV 2 , 



Af =(0.54 ±0.21) • 10~~ 3 GeV 



(56) 



We choose a g = a 3 /2 1 ^ 6 = 0.65 GeV -1 and a s = 0.5 at 
the scale 1 GeV which results in the following probabil- 
ities for the quark-gluon components within the nucleon 
state IH 



35 
105 



P 9± =i^ m NP^i) 2 - 0-15 

(Af) 2 + (Af) 2 



16g 2 



0.185. 



(57) 



RESULTS AND DISCUSSION 



Now that we have models for the nucleon LCWFs, it 
is straightforward to evaluate the matrix elements of the 
four-fermion operators Q± and constrain the momen- 
tum fraction dependence of the corresponding higher- 
twist correlator Q^ p \ The distributions Q±(£) defined 
by Eq. (I32p possess the following support properties 



C±(0 =0(-&)0(-6)0(6W&) 

x d(i - e 2 - 4 )<z ± (-6 , 6, -6 , 6) ■ (58) 



.01 
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-0 . 02 
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. 0005 




-0 . 0010 



-0.0015 



FIG. 2: The nucleon twist-four distribution Q^ p \x) multi- 
plied by a% (solid line). The dashed and dotted lines show 
the contribution of three-quark and and quark-gluon wave 
functions, respectively. The lower panel is a blow-up of the 
high x region. 



Here the functions q±(£) are expressed in terms of in- 
tegrals involving the nucleon wave functions, see Ap- 
pendix [B] for explicit formulas, while below we quote ex- 
pressions which correspond to the ansatzes (|49l) and (|55l) . 
The structure of the Fock expansion corresponds to the 
decomposition of the twist-four distributions q± into the 
following three components 



?±(0=9±'(0+3± + (0 + 



(59) 



Each term in this sum corresponds to the contribution 
of the pertinent multi-parton component of the nucleon 
wave functions, i.e., three-quark and quark-gluon, respec- 
tively. Making use of the results derived in the previous 
section, one finds the following explicit momentum frac- 
tion dependence for the distributions q± , 



Q- q (0 =c 3q xi{0 [(4 - 3(6 + £ 4 )) 2 + (5 - 36) (5 - 36) 
«+*(0=ca 9 xi(0 (1 + 3^0(1 + 36), (60) 

where 

1-6-6 



xi (0 =6666 



6 + 6 



(61) 
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and the overall normalization constant being 

560 

C3q = r-Sq 



87ir 2 al ■ 



(62) 



For the four-parton quark-gluon functions q± , q± one 
gets 



where 



Xa(0=Xi(0(l-&-e4) 3 



(63) 



(64) 



and 



-P, 



560 



\ 2 
A 3 



700 



n~ = — P 

°9t r 3t „2„2 



+ _p 280 



3 A| + A3 
A 3 



^2 + ^3 



6 A 4 A 
5 A 2 - -A 3 



T On 



while = 0. 



(65) 



Furthermore, making use of Eq. (|42p one obtains after 
some algebra the following representation for the function 
Q p (x), a; > 0: 



Q P (x) 



1 



1 



1—X 

2 1 dtij]ag(x/Q$+(x,{,Z,x)+ 



1 
2 



1 r x - x -tdr) 



V 



x + ( 



\q+far),t,x + £- 77) - ~q + (x,Z,£,x)) - q-(x,r),£,x + £- 77) 

r) x + £ 



?7 

77 + x_ 



T) + X £ 



q-(ri + x,ri,£,£ + x) 



(66) 



where 

$±(0 =^(1 + ^14^23) (l + Pi3^ 24 ) q±(0- (67) 

Performing the final integration is straightforward and 
one can obtain a closed analytical form of Q P (x) (how- 
ever, the resulting expression is quite long and in order to 
save space it will not be displayed here). The twist-four 
distribution is displayed in the upper panel of Fig.[2J The 
dashed and dotted lines correspond to its three-quark and 
quark-gluon components, respectively. Both of them ex- 
hibit a global minimum at x ~ 0.4. In the lower panel 
of Fig. [21 we blow up its high- a; region to demonstrate 
the node structure of the three-quark contribution. As 
x — > 1 the four-parton quark-gluon component of Q p (x) 
is suppressed by the decay factor (1 — x) 3 with respect to 
the three-quark component. At the same time the twist- 
four distribution Q p (x) is enhanced in comparison with 
the twist-two parton densities calculated within the same 
model, Qp(x)/u p (x) ~ log(l — x) for x — > 1. 

Our predictions for the twist-four correction i?* w ~ 4 to 
the Cahn- Gilman formula is shown in Fig. [3j In order to 
make an comparison with the results of Ref. |14| easier, 
we display ijf -4 for Q 2 = 4, 6, 8, 10, 12 GeV 2 . It turns 
out that our prediction for i?^ w ~ 4 is roughly twice as 
large as that of Ref. [Ti[ with the minimum of the func- 
tion being slightly shifted towards lower x' (i.e., from 
x' ~ 0.7 to x' ~ 0.6). Note that the x-dependence of the 
twist-four contribution is much better determined than 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

x — 2x 



FIG. 3: The estimate 7i' w_4 as a function of the Bjorken x for 
different values of Q 2 . The curves from the bottom to top cor- 
respond to the values Q 2 — 4, 6, 8, 10, 12 GeV 2 , respectively. 
The experimental accuracy of SoLID is ±0.005 for i?J™ -4 at 
an average Q 2 of 3.3 GeV 2 and (x) = 0.34. 



its normalization: The three-quark component of the nu- 
cleon wave functions is constrained by the existing exper- 
imental data (parton densities and nucleon form factor, 
[HI), but the ansatz ([55)1 for the quark-gluon wave func- 
tions has to be regarded as an exploratory estimate (see 
Ref. for a discussion). Nevertheless, since for large 
x 1 the contribution due to the quark-gluon components 
of the wave functions are strongly suppressed, see Fig. [21 
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we believe that for x' > 0.7 our estimate for R\ w ~ i (x') 
should be rather accurate. That is, the function Ri(x') 
has to change sign around x' ~ 0.8. We also checked 
that our result, once we compute its Mellin moments, 
are in good agreement with earlier calculations of higher 
twist corrections to the first moments of structure func- 
tions (3 Hi. 



VI. CONCLUSION 

Parity-violating deep inelastic scattering is a process 
of fundamental importance and, therefore, will be inves- 
tigated by ever more precise experiments. It is sensi- 
tive to physics beyond the Standard Model as well as to 
specific aspects of strong interaction dynamics, encoded 
in higher-twist correlators. To disentangle both, the IE- 
dependence of the twist-four contribution must be known 
precisely which seems to be in reach with present day 
techniques. The task of determining these higher-twist 
contributions has a certain urgency in view of upcoming 
JLab experiment SoLID [To| . In the current study we cal- 
culated the twist-four correction to the leading contribu- 
tion di to the parity violating asymmetry by determining 
matrix elements of light-cone four-quark operators fl6| . 
We found that within the framework of light-cone wave 
functions, the estimate for twist-four correlation func- 
tions has similar features as found in a recent calculation 
within the MIT bag model [lij]. The size of the correc- 
tion Ri is about twice as large in our calculation and the 
form differs slightly, but these differences might well re- 
flect the present day theoretical uncertainties of such cal- 
culations. The size of the twist-four correction we obtain 
is borderline. It has to be taken into account to improve 
the sensitivity of SoLID for New Physics, but it does not 
seem to be large enough for SoLID to test our predic- 
tion. However, as mapping out the running of sin 2 9\y is 
one of the fundamentally important experiments we are 
optimistic that still more precise experiments will be per- 
formed in future, which should then be sensitive enough 
to observe the higher-twist contributions we analysed. 
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Appendix A: Light-Cone expansion 

In this Appendix in order to make the paper self- 
consistent, we spell out our notations and conventions 



that we used to perform calculations of hadronic matrix 
elements in the body of the paper. 

For an arbitrary four- vector a M we define the light-cone 
coordinates as 



1 



a+ = -^=(a° + a 3 ) , a- = -j=(a° - a 3 ) . 



a = a 1 + ia 2 . 



a = a — ta~ 



(Al) 



We find it convenient to pass from four-dimensional vec- 
tors to two-dimensional matrix notations for all tensors. 
For a vector we introduce the matrix a — a^a^, where 
<r" = (!,<?), 



n _ ( V2a- —a 



(A.2) 



In the Weil representation the Dirac 7— matrices has the 
form 



I 



4 0' ^ V -o l 



cr' 



7 



-I 
I 



with 7 5 = i7°7 1 7 2 7 3 . In the two-component notation 
the Dirac spinors read 

q = fe) , g = g t 7° = (^ffr), (A.3) 

where gfm = — (1 ± 7s)<Z are components with posi- 
tive/negative helicity, respectively. The two independent 
light-like vectors 



1 



V2 



(1,0,0,-1), n» 



1 



V2 



(1,0,0,1), (A.4) 



n 2 = n 2 = 0, n ■ n — 1 can be parameterized in terms of 
two auxiliary Weil spinors: 

which read explicitly 

k = 2 1 / 4 , Ha = 2 1 / 4 (; 

A^V^- 1 ), ^=2V4^. (A.6) 

The following rules allow to raise and lower spinor indices 

A" = e a/3 X/3 , \ a — A^e^Q , A" = X^e 0a , Aq = e d/ jA^, 

with the antisymmetric Levi-Civita tensor having only 
the following nonzero components 

e 12 = e 12 = - £i 2 = -e i2 = l. 
The auxiliary spinors A and /i are normalized as 

{flX) = fi&X & = -(^p) = +V2 (A.7) 
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and are used to project out the "plus" and "minus" com- 
ponents of the fields. For fermions, we define 



X+ = XdA Q , 



(A., 



In the same fashion the light-cone decomposition of a 
vector (e.g., gluon) field takes the form 



A_ X a Xa + A + 

H'CtH'dt 



V2 



A 
V2 



fJ>a Xa 



The "plus" spinor fields ip+ , \+ and transverse gluon 
fields A, A are assumed to be the dynamical degrees of 
freedom in the light-cone quantization framework. While 
the "minus" fields "0-, X—> A- can be expressed in terms 
of these with the help of equations of motion. Finally, 
we use the gauge A + = 0. 

The good components of the quark field have the fol- 
lowing canonical expansion 



dp+ d 2 p 



v/2pT (2tt) 3 
dp+ d 2 p ± 
v/2pT(2^) 



- ip - x h(p) + e+^ x d\{p) 

- tp - x b t (p) + e +ip - x d\(p) 
(A.9) 



in terms of the annihilation operators of quark and an- 
tiquark of positive (negative) helicity b^m , dfm , respec- 
tively. They obey the standard anticommutation rela- 
tions 



{b x (p),bl,(p')} = {d x (p),dl,(p')} = 
= 2p + {2ir)H x . x ,5(p + -p' + )5^( Pl 



Pa 



(A.10) 



Similarly the expansion for the dynamical transversely 
polarized gluon fields A and A reads 



7/ % rx f dk+ d 2 k± , . 

, , n rx f dk + d 2 k^ , 
A ^ = ^M(2^^ 



- i k 



x at(k) + e +ik - x a\(k) 



e~ ik - x ai(k) + e +lk - x a\(k) 

(A.ll) 



Here and below A = ^ a t a A a are matrices in the fun- 
damental representation of SU(3) and t a are the usual 

generators, normalized as ti(t a t b ) = ^ ah ■ The creation 

and annihilation operators obey the commutation rela- 
tion 



a b x (p),(a b ;,(p')y 



(A.12) 



= 2 P+ {2nf8 x , xl 8 bb '5{p + -p' + )S^(p ± -p' ± ) . 

As mentioned above, bad (i.e., "minus") components can 
be expressed in terms of the dynamical fields using QCD 
equations of motion. 

Appendix B: Distributions q± 



As discussed in the main text in Sect. |Vj we represent 
the twist-four distributions q±(£) as shown in Eq. (|59p . 
We remind here that the arguments £ = (6,6, 6, £4) are 
subject to the constraints, < 6 < 1 and 6+6 = 6+6- 
A straightforward calculation of its components <z± 9 (0, 
q±((i), 9± T (6, arising from three- and four-parton Fock 
states of the nuclcon, yields the following expressions in 
terms of the LCWFs introduced in the main text, 



= 

d q (0 = 



(™3/Af) 2 



9 (£2 +&)(!- 6- £4) 
4 {na 3 f N ) 2 



0(6 , 1 - 6 - 6, 6)0(6, 1 - 6 - 6, 6) , 
1-6-6,6,6)0(1-6-6,6,6) 



9(6 + 6)(i-6-6) 
+ (0(6, 6, 1 - 6 - 6) + 0(1 - 6 - 6, 6, 6)) (0(6, 6, 1-6-6 



(i -6 -6,6,6: 



(B.l) 



where is given by Eq. (|48|). Next, we got that = and 



32(87r 2 a2) 2 f 1 dx 2 dx 4 



, r^ 1 ~ ^2 -6 - x 2 - x 4 ) < S (6, 3:2,6,^4) 4> g {x2, 6, 6,2:4) + -0 s (6,a;2, 6, ^4) 

0(6 + 6) JO x 2 x 4 



+ <f>g{%2, 6:6:^4) 05(6,^2,6,^4) - 20 S (X 2 ,6,6,2 ; 4) 



(B.2) 
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Finally, 

tf(0 = 



9^(0- 



oL i e\ / — ~T :,5 ( 1 -C2-'C4-a;2-a;4)'|i/'g 1) (6,a;2,6,a ; 4) (fa, &2, £4, x 4 ) + 5^ 2) (£ 2 , £4, x 2 , au) 

"Hs2 + S4j Jo ^2 £4 I L 

+ (£l ) 6 , X 2 , X4 ) (6 , £4 , £2 , X 4 ) + 5lpW (£2 , X 2 , £4 , X 4 ) 

32(8A 2 ,) 2 f 1 dx 2 dx 4xn A r ^ 1 / * a 

3(6 + £4) Jo x 2 xi y 1 4 

+ ^9(^2, £l,6, ^4) </>g(6, 3C2, £4, 2:4) - 2^(^2,6) £4, Xi) - ^ 1) (x 2 ,^2,U,X 4 ) ~ j1pl 2) (x2,U,^2,X4) 



X lp( 2} (x 2 ,&,£l,X4) ~ 1p^ ) (x2,^1^3,X 4 ) 1p (2) (x 2 ,U,(,2,X4 : ) + 2tp^(x 2 , 6, ^4, Xi) 



(B.3) 



where 



tp g [x 1 ,x 2 ,x 3 ,X4,) = ^ 1) (x 1 ,x 2 ,x 3 ,x 4 ) - ip t X 2) (x 3 ,x 1 ,X2,x 4 ) . 
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